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We discuss heavy quark production in High Parton Density QCD in quasi-classical approximation 
and including \ow-x quantum evolution. We also consider an alternative approach based on the effect 
of pair production in external fields. 
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O ' I. QUASI-CLASSICAL APPROXIMATION 

<n : 

Production of a quark q and antiqiuark q pair at high energies is characterized by two time scales: production 
time Tp and interaction time r m t. In pA collisions in the center-of-mass frame of the qq pair the production time is 
rp ~ l/(2m), where m is a quark's mass. In the nucleus rest frame, this time is Lorentz time-dilated by E g / (2m) where 
Eg is energy of a gluon in a proton from which the qq originates. In terms of the Bjorken variable X2 = (nriT I \/s)e~ y 
the production time is Tp ~ 1/ (2Mx2), where M is a nucleon mass. At RHIC this corresponds to the production 
time Tp ~ I5e y fm. On the other hand, the interaction time is Tint — Ra — 7 fm. Clearly, in the limit y 3> 1 one can 
H ' consider the pA interaction as an instantaneous process. This implies in particular, that in the configuration space 
^ | the transverse coordinates of quarks and gluons do not change over time of interaction in which case color dipoles 
ON • diagonalize the interaction matrix. 

The time structure of heavy quark production discussed in the previous paragraph can be taken into account in 
the light-cone perturbation theory. It is convenient to use the nucleus rest frame. In Ref. a model was considered 
in which qq pair appears as a result of the following chain of fluctuations in a fast proton: q v — > g — ► qq, where q v is 
a valence quark and g is a gluon. Thus, the heavy quark production amplitude consists of contributions from three 
different time orderings: (i) valence quark interacts with the nucleus and then emits a gluon, (ii) gluon interacts with 
the nucleus and then emits a qq pair, (iii) qq interacts with the nucleus. Correspondingly, there are six different time 
orderings in the cross section [y • The light-cone wave function for each time ordering was calculated in Ref. Q . In 
addition, one also has to calculate the propagator of a parton system for each time ordering. The propagator sums 
Oh. up multiple scattering of a dipole in a nucleus. In the quasi-classical approximation this amounts to resuming higher 
order corrections in a s enhanced bvpowers of large nuclear length ctgA 1 / 3 ~ 1. Explicit expressions for all these 
propagators can be found in Ref. 0, Q • 

The final result for the double-inclusive quark-anti-quark production cross section in pA collisions in the quasi- 
| classical approximation is Q 
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where k l , k 2 are the quark and anti-quark transverse momenta, y is quark's rapidity, a is a fraction of gluon's energy 
carried by a quark, b is impact parameter. $^ and Sy — 1,2,3) are light-cone wave functions and propagators 
for different time orderings in configuration space. Integrations in Q go over all transverse coordinates in a proton 
wave function. 

The complicated result of Ref. can be simplified for the aligned jet configuration: when quark or anti-quark 
takes away most of the gluon's energy (a -C 1) I]- In that case the gluon production time becomes much larger than 
that of qq which allows factorization of the light-cone wave functions into those of gluon <& g ^ q q and quark-anti-quark 
<& qv —>q v g. Then the differential cross section for single inclusive quark production in pA collisions reads |3j,|^,|5| 
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where I assumed for simplicity that the dominant contribution comes from interaction of the qq pair with the target, 
while rescatterings of the gluon and the valence quark are neglected. However in general, they must be taken into 
account as well. The ligh-cone wave functions in (0 are given by 

, a s C F (ax + (1 - ops') ■ {ay + (1 - a)x') 

— 7r^ (ax + (1 — a)x') z (ay + (1 — a)x') z 



$g^qq{z,X,X_',a) = — TV 2 
7T 



(x — x 1 ) ■ (y — x') 
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x[a 2 + (l-a) 2 }+ K (\x-x'\m)K (\y-x'\m)j , (4) 

where Q s is the saturation scale, m is a quark mass. 
The gluon saturation scale Q s is given by 

Q* = ^^pT(b)xG N (x,l/x*), (5) 
Of 

where T(b) is a nucleus profile function and p is the nuclear density. The saturation scale measures the strength of 
the nuclear color field at high energies F ~ Ql/g (see reviews jy, Q)- 



II. INCLUDING QUANTUM EVOLUTION 



When the interaction energy becomes such that a s y ~ 1 it becomes necessary to include the effects of quantum 
evolution at small x into the qq pair production cross section 0, Q ■ This is done by first generalizing the quasi- 
classical model to include possible gluon emission from valence anti-quark. Thus, we get a formula for the case of 
qq production in dipole-nucleus scattering. Strictly speaking, our results would then only be applicable to particle 
production in deep inelastic scattering. However, our results below may still serve as a good approximation for gluon 
production in pA collisions . If the transverse coordinates of the quark and anti-quark in the incoming dipole are 
denoted by z and z_ x correspondingly with z 01 = z — z_ x , we write 

d^hd^dydad^- 01 ' - 4(27r)6 J d yi d y 2 e 
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with generalized propagators 0- 

The inclusion of quantum corrections due to leading logarithmic (resumming powers of a s y) approximation in the 
large- N c limit is done along the lines of (see also for a review) using Mueller's dipole model formalism Since 
the integration over rapidity interval separating the quark and the anti-quark in the pair does not generate a factor 
of the total rapidity interval Y of the collision (i.e., does not give a leading logarithm of energy), the prescription 
for inclusion of quantum evolution is identical to the single gluon production case. We first define the quantity 
n i(ikiT z_i] Zjo> i z.ii]Y — y), which has the meaning of the number of dipoles with transverse coordinates ZQ>,Zy at 
rapidity y generated by the evolution from the original dipole Zq,z_i having rapidity Y. It obeys the dipole equivalent 
of the BFKL evolution equation 0, 0] 

ni(z ,z 2 ;z ,,z lt ;y) + ni(z 2 , z x ; z , , z v ;y) 
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with the initial condition 

ni(z o ,z 1 ;z o ,,z 1 ,\y = 0) = S(z - Zq,) 5^ - z v ). (8) 

If the target nucleus has rapidity 0, the incoming dipolc has rapidity Y, and the produced quarks have rapidity y, the 
inclusion of small- a; evolution in the rapidity interval Y — y is accomplished by replacing the cross section from © 

byS 

(z 01 ) - J d 2 z , d"z v n^z^^Y - y) ——^— {gvv ). (9) 



d 2 k\ d 2 k 2 dy da d 2 b 



Inclusion of evolution in the interval between and y is accomplished by replacing the Mueller-Glauber rescattering 
exponents according to the following rule 

e-ifeo-SaJ'OjMVIaa-^IA) ^ Nix^x^Y) (10) 

where N(xq,x_i,Y) is the forward amplitude for a quark dipole x Ql x l scattering on a target with rapidity interval 
Y between the dipole and the target. (We refer the interested reader to the Ref. Q for explicit expressions of the 
resulting propagators.) The forward scatetring amplitude obeys the following evolution equation Q 
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with the initial condition 

N(x ,x 1 ,Y = Q) = l- e -3fe>- £ i )2Q ' ln (Vli£o-£il A ). (13) 
The final result for the double inclusive qq production cross section including small- a; evolution effects reads 

do 1 f 
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x Y (" 1 ) fe+ ' $ yfei~^fe>£2-2fe;y 1 -^ J i2 _ -i ;a ) E ij(z.i'X2,Zk',y v y 2 ,zi', a ,y)- ( 14 ) 

ij = l fc,2=0 

Integrating over one of the quarks' transverse momenta we obtain the single inclusive quark production cross section 
"Uoi) = m [ d 2 z Ql d 2 z v ni{z^,z 1 ]z Wl z v ]Y -y)d 2 xid 2 x 2 d 2 yie~ l - { - 1 ^i ) 
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da Y Y ( _1 ) fe+!$ yfei -l k ,x 2 -z k -,y_ x -z,,z 2 -z t ;a) ^{x^x^z^^x^^-a.y). (15) 
~' i,j=i fc,;=o 



These equations generalize the results Refs. 0,0, H,0 
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III. SCHWINGER MECHANISM 



Alternative approach to particle production in high energy hadron and nuclei collisions was advocated in Refs. [l3l 
llH . It was argued there that the particle production from the Color Glass Condensate can be considered as a vacuum 
instability in strong background longitudinal color field. It has been known for a long time that strong longitudinal 
color fields may be important for particle production in hadronic reactions |l5j and in particular, in heavy-ion collisions 
|IcL ITU ITU H^. I20I l2lj . However, the structure of these fields remained unknown due to their non-perturbative nature 
characterized by large coupling a s ~ 1. The Color Glass Condensate eliminates this problem by introducing the hard 
scale Q s 3> Aqcd so that the classical color fields can be in principle calculated as a series in small coupling and high 
density of charges a s {Q 2 s ) p ~ 1, where p oc A 1 / 3 . Existence and dominance of longitudinal fields in color glass was 
realized in [T3L FLU ] and later in ■ 

Moreover, according to the parton model, particles in a parton cascade are strongly ordered in energy. This implies 
that the production times of a partons in a cascade are also strongly ordered. Therefore, it makes sense to consider 
motion of a slow parton in an almost static field of faster partons. One can also argue that the transverse sizes of 
partons are strongly ordered as well, at least in the leading logarithmic approximation, which makes the longitudinal 
chromoelectric field E the dominant in the target rest frame |14| . 

Consider for simplicity motion of a scalar charged particle <fi in the background field 

(d M - igArf = 0. (16) 

The field is given by E z = —d-A+ = const. Although A+ is a function of only x- Eq. (|16J) cannot be solved by 
separation of variables since the initial condition depends on both x+ and X-. Indeed, a fast hadron interacts with a 
target located at, say, z = t = 0. Instantaneously, the parton cascade looses its coherence. Therefore, we have to solve 
equation of motion i|16[l with the constraint that the potential A + vanishes at z = 0, i.e. explicitly depends on both 
lightcone coordinates x+ and X-. Thus, we are looking for the solution in the form 4>{x) = e~ ipx ' XJ - . Working in 
the WKB approximation \d+Sd-S\ ^> \d+d-S\ we reduce ijTTfl) to 

-2d+S(d-S - gA + (x-)) + p\ + 2gaE z = 0, x+ > x_ , (17) 

where <9+ = and cL = Eq. (|16|l is a Hamilton- Jacobi equation of motion of a charged particle in the 

background field A + = A- (x + )9(x + — xJ). We can write and solve the Hamilton equations of motion for the 
canonical momenta p_ = — d-S, p+ = —d+S as a functions of the light-cone coordinates. The result is [Tj] 

P+ = -gA + {x-) + gA+{x+) + P ° + , (18) 
dx+ 

(p\+gA + {x + )f 

gA+{x - )dX - + J dX+ P l + gi + (x + ) - ^ 

Eq. I|19|l coincides with the equation of motion of a classical test particle of mass p± in the external field A+ {x+ ) . In 
other words, the test particle effectively moves under the action of the longitudinal electric field E z = — A',(x+). 

Eq. (|20|l gives the action of the test particle along the trajectory l|19l) . Its imaginary part arises from the pole 
in the integrand of the second term in the right-hand-side of (|20() . Integration around the pole in the plain of 
complex x + yields the imaginary part. It can be calculated replacing the denominator in the first integral in l|2(J|l by 
Im(p° +5^4+) _1 = ±(iir /2)8{p a + + gA + ) according to the Landau-Cutkosky cutting rule. Additional factor of 1/2 
arises due to the condition x + > x_ . Define 

r = x + u;, A + (t) = -^/(r), 7=4^' ( 21 ) 

uj gE 

where w is a typical frequency of the external field and E(t = 0) — Eq. With this definitions we obtain 

l m5 = Im / 4-^-- = £r „,Zr ^ ■ ( 22 ) 



*- = 4/ 73 , dX +„ (19) 



gE 01 -f(r) 2 3 S /'(/- 1 (7)) 

The imaginary part of the action l|22|l corresponds to the pair production. 

The physical meaning of the adiabaticity parameter 7 introduced in (|21|l is clear: 7 = for the static field, while 
7 ^> 1 for rapidly oscillating one. Since gE^ ~ kf , and to — ki t -, with subscript i denoting a parton making up the 
background field, we have the following estimate 

K i,+ 
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This estimate implies that due to the strong ordering of the light-cone momenta in the partonic cascade at r = 0, 
the emission of the gluons is determined by small values of 7 or, in other words, by constant electric fields, in which 
A + (x + ) = E x + . Using this formula in yields the result for pair production probability [23| 

" eXP {" gE J' (24) 

We would emphasize that (|24[1 is not analytic at small coupling which is a manifestation of its non-perturbative nature. 
At small transverse momenta it has a finite limit even at m = 0. Therefore, the proposed method of calculating the 
pair production at high energies allows access to the kinematical region of small p± where the perturbation theory 
breaks down. 

As later times the background field decays due to the effect of screening by produced pairs which results in modi- 
fication of the spectrum [l3l lh"H . 



IV. PHENOMENOLOGICAL APPLICATIONS 



Equations l|14|) , l|15() and (|24() have important phenomenological applications for studying the dense partonic system 
in p(d)A and eA collisions. Observation of hadron suppression in the nuclear modification factor measured in dA 
collisions at forward rapidities at Relativistic Heavy Ion Collider (RHIC) |2^] signals the onset of the nonlinear 
evolution of the scattering amplitude for light hadrons [2J l2Sl |29| . Due to a large mass, the impact of nonlinear 
evolution effects on the heavy quark production is shifted to higher energy and/or rapidity. It was estimated in j^lj 
using the fcr-factorization approach (see also |31|) that one can expect a significant deviation of the open charm 
production cross section from the perturbative behavior already at pseudo-rapidity rj ~ 2 at RHIC. Due to the heavy 
quark production threshold one expects that the total multiplicity of open charm scales as iV C oii at lower energy and /or 
rapidity whereas at higher energies and/or rapidities the scaling law should coincide with that for lighter hadrons |30| . 
i. e. open charm multiplicity should scale as N part [13 due to high parton density effects. Therefore, to be able to 
compare predictions of CGC with the data reported by RHIC experiments and to make predictions for the possible 
upcoming pA run at the Large Hadron Collider (LHC), it is important to perform a calculation of an open charm 
production within the more general approach developed in this paper. Our final results l|14fl. (|I5|I and (|24|l allow one 
to describe open charm transverse momentum spectra at different rapidities and center-of-mass energies, allowing for 
a complete description of RHIC and LHC data. Since the saturation scale Q s is expected to be even higher at LHC 
than it was at RHIC, the CGC effects on heavy quark production at LHC should be even more significant. 
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